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Abstract. Using active gel theory we study theoretically the properties of the cortical actin layer of an- 
imal cells. The cortical layer is described as a non-equilibrium wetting film on the cell membrane. The 
actin density is approximately constant in the layer and jumps to zero at its edge. The layer thickness is 
determined by the ratio of the polymerization velocity and the depolymerization rate of actin. 
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1 Introduction ing with a large number of proteins among which myosin 

I molecular motors play a major role. Myosin motors, as- 



■ Living cells maintain and change shape, adhere, spread, 
^ ' divide and erawl through the agency of a dynamic, fila- 



sembled in minifilaments, consume free energy through 
the hydrolysis of ATP molecules and can produce work. 



■ mentous scaffold known as the cytoskeleton IT]. Although „i-t i ■ m ■ ■ -m 

$H ■ L5y binding to the actm nlaments, myosin mmitilaments 

this structure is made up of many different constituents ., , ... ,. , 

create contractile stresses m the actm gel. 

the task of stress generation, with crucial consequences 

for the mechanical properties of cells, lies primarily with The crosslinked polymer network formed by actin and 

the acto-myosin component 0. This substructure consists its associated proteins differs profoundly from more famil- 

of a meshwork of semi-flexible actin filaments interact- iar thermal-equilibrium physical or chemical gels because 



^ on leave from Department of Physics, Indian Institute of- of the sustained energy dissipation by the molecular mo- 
Science, Bangalore 560 012 tors and some other proteins. It is usefully viewed as a 
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state of active matter [3] known as an active polar gel [H 
[51[SJ[71[S] - "active" referring to the steady consumption 
of free energy at the scale of individual components, i.e., 
the actin-bound molecular motors, and "polar" to the ori- 
entable, directed character of the actin filaments. In par- 
ticular, each filament carries a distinction between its two 
extremities, as does each constituent monomer. We will 
therefore refer to the plus and the minus end of a fila- 
ment. 

The filaments of the cytoskeleton undergo constant as- 
sembly and disassembly: Each actin filament polymerizes 
preferentially at the plus end and depolymerizes at the 
minus end. This process is called treadmilling, and is reg- 
ulated by a multitude of accessory proteins. Through de- 
polymerization, filaments disintegrate, and through poly- 
merization new filaments are generated, maintaining a 
nonzero gel mass on average. Generation of new filaments 
is assisted by nucleating proteins, such as formins that act 
as seeds for actin polymerization and then stay attached 
for some time to actin plus-ends, where they promote the 
addition of monomers. The Arp2/3 complex instead binds 
to existing filaments and remains attached to the minus 
end of newly created actin filaments. 

In this work we study consequences of the interplay 
of actin polymerization and active contraction of actin 
gels through a hydrodynamic description, that captures 
the generic behavior of materials on large length and time 
scales. Such a description starts with the correct choice 
of slow variables [9] relying broadly on conservation laws, 
broken continuous symmetries, and order-parameter modes 
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near a continuous phase transition. This approach has 
been successfully extended to active systems [51B1 ITU1ITT] . 
held in stationary states far from thermal equilibrium by 
the sustained dissipation of free energy. Despite the very 
broad range of length scales, from microns to kilometers, 
on which organized active matter is seen, there is a degree 
of universality in its hydrodynamic properties, classified 
by the type of broken symmetry and applicable conser- 
vation laws: examples include waves without conventional 
inertia, anomalously large number fluctuations, and a ten- 
dency of instability of quiescent states towards sponta- 
neous flow [3]. 

We focus on the fact that a substantial fraction of 
the actomyosin is located in a layer, known as the actin 
cortex [T21II3], adjacent to the plasma membrane of the 
cell, see Fig. [T^. This localization of one component to 
the vicinity of a wall is reminiscent of the physical phe- 
nomenon of wetting. However, prevalent explanations of 
the cortical actin profile do not take advantage of this 
analogy, relying instead on an imposed spatial separation 
between the zones of polymerization and depolymeriza- 
tion. In this paper we show that the active contractil- 
ity alluded to in an earlier paragraph can drive a tran- 
sition to a steady state maintained by a polymerization- 
depolymerization process, sec Fig. [TJ), in which the poly- 
mer concentration has a profile very similar to that of a 
wetting layer [14j . The cortex thickness and actin density 
profile emerge naturally from our treatment. Recall that 
in thermal equilibrium systems wetting arises through the 
selective attraction of a component to a surface. In the fol- 
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I"").' 2 Hydrodynamic description 

L 



(»• Wc now present the hydrodynamic equations governing 
'-'k; the dynamics of an assembling active actin gel. We assume 
that the gel is assembled by actin polymerization at a sur- 
face. In a cell or in an in vitro experiment, the polymeriza- 
tion is promoted by nucleating proteins, such as proteins 
of the formin family, see Fig. [Ij). Away from the surface, 
the actin gel assembles by elongation of existing filaments 
or by nucleation of new filaments. The gel disassembles be- 
cause of monomer removal at filament minus-ends or by 
severing of the gel filaments that produces small filaments, 
which diffuse into the solution. Myosin molecular motors 
assemble into small filaments that act as cross-links, which 



Fig. 1. a) Confocal fluorescence microscopy image of a HeLa 
cell with F-actin labeled by life-act ruby. Below the plasma 
membrane a well-defined zone of high fluorescence represents 
the actin cortex. Image with courtesy from M. Fritzsche and 
G. Charras. b) Schematic illustration of the cortex dynamics 
considered in this work. Nucleation promoting factors located 
at the plasma membrane nucleate new actin filaments or assist 
elongation of existing filaments at rate k*. This generates a 
flux of polymerized actin Vppo, where po is the gel density at 
the surface and Vp = is the polymerization velocity, 5 being 
the size of an actin monomer. Filaments disassemble at a rate actively generate mechanical stress in the filament net- 



fed anywhere in the cortex. Not shown is the growth of actin work. The network is permeated by a solvent containing, 

filaments away from the membrane. Molecular motors act as in particular, unbound motors and actin monomers. We 

active cross-links and generate active stresses in the cortex. limit ourselves to the case where the exchange of motors 

between the actin network and the solvent is so fast that 

lowing we show that contractility, although quite different ,i , i -ru * j j u -ti 

° 6 H we can assume them to be equilibrated and where the mo- 

from an attractive potential, plays a similar role in driving ^^^^^ diffusion is fast enough that the concentration of free 

a condensation at the plasma membrane. To obtain this ^^^^^^^.^ ^^j^^^^ constant. In a first approximation, 

result we generalize the hydrodynamic description of ac- ^^^^ concentration of motors bound to actin is then propor- 

tive gels to include crucial nonlinear density dependences, ^-^^^j ^j^^ j^^^l ^^^-^ concentration. We have verified 

The rest of this paper is organized as follows. In sec- numerically that our results do not change significantly 

tion II we present the hydrodynamic equations for poly- when motor binding and unbinding and motor diffusion 

merizing active gels. In Section III, we study an active gel are taken explicitly into account, at least in the case when 

polymerizing at a surface, which yields a description of the the gel does not impede the diffusion of motors, 

actin cortex of animal cells. The properties of the gel can As mentioned in the introduction, actin filaments are 

be understood in terms of active wetting and dewetting polar objects. Consequently, an actin gel can present a 

by filament assembly and disassembly. macroscopic vectorial order represented by a polarization 
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field p. In the present treatment we ignore the subtleties 
associated with actin polarization and consider in the fol- 
lowing the case of an isotropic gel. Our detailed calcula- 
tions will be for the case where the actin concentration 
varies only along one direction, where transverse varia- 
tions in p play no role. 

Actin polymerizes on a planar surface located in the 
plane z = 0. The gel and the solvent are restricted to 
the half-space z > 0. If the gel is homogeneous in the di- 
rections parallel to the surface, its properties and in par- 
ticular its density only depend on the z coordinate. We 
will explore elsewhere the dynamics of variations in the 
xy-plane. 

2.1 Conservation laws 

The hydrodynamic description is based on conservation 
laws for mass and momentum. Mass conservation of the 
gel and the solvent read 

dtPg + daPgVg^a = -kdPg 

dtPs + daPsVs,a = kdPg ■ (1) 

Greek indices denote the three spatial directions x, y, and 
z and we have adopted Einstein's summation convention. 

In the above expressions, pg and ps, respectively, de- 
note the densities of the gel and the solvent, while Vg and 
Vg are the corresponding velocities. Degradation of the gel 
due to depolymerization and severing occurs at a constant 
rate kd- The gel is produced on the one hand by growth 
of existing filaments, and on the other hand by the nucle- 
ation of new actin filaments. There is no significant spon- 
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taneous nucleation of new filaments at concentrations of 
monomeric actin present in cells. Instead nucleation pro- 
moting factors regulate the generation of new filaments 
and their elongation. Members of the formin family and 
the Arp2/3 complex are important examples of nucleation 
promoting factors P31[TC1[T7| . Formin is located directly 
beneath the plasma membrane, where it generates new 
filaments. The Arp2/3 complex needs to bind to existing 
filaments before it can act as the seed of a new filament. 
Still, it is predominantly localized in the vicinity of the cell 
membrane. We therefore consider here for simplicity that 
the actin gel growth occurs only at the surface. Therefore, 
the filament density far away from the surface vanishes. 
We account for elongation and nucleation on the surface 
by a boundary condition on the gel fiux at z = 0. Ex- 
plicitly, we write PgVg,z\z^Q = VpPo, where Vp is the poly- 
merization velocity and po the gel density at the surface, 
which is imposed by the density of nucleation promoting 
factors. 

2.2 Constitutive relations 

To fully specify the behavior of the actin gel, we must 
provide constitutive equations. They link the generalized 
thermodynamic fluxes to the generalized thermodynamic 
forces. We follow here closely the approach of Callan- Jones 
and Jiilicher [TS] to describe active permeating gels. The 
fluxes are in our case the gel velocity, the time derivative 
of the strain in the gel, the relative current between actin 
gel and solvent 3 = Pg (vg — v) , where v is the center of 
mass velocity, the deviatoric stress tensor a and the rate of 
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ATP consumption. The generalized forces are the gradient 
in the relative chemical potential p, of the gel and the sol- 
vent, the center of mass velocity gradient Vv, the partial 
stress of the gel crS and the activity of the system. The rel- 
ative chemical potential is given by /I = fig/nig — /ig/mg, 
where fig and /is are the respective chemical potentials of 
the gel and the solvent and mg and the masses of the 
solvent and actin monomer molecules. Active processes are 
eventually driven by the hydrolysis of ATP into ADP and 
inorganic phosphate P^, with chemical potentials (iatp, 
Madp, and /ip, respectively. We express the system's ac- 
tivity through the difference Afi — /iatp ~ A'adp ~ Mp- In 
the following Afi is considered as constant in space and 
time. We do not consider any further the corresponding 
thermodynamic flux which provides the rate r of ATP 
consumption in the system. 

The constitutive equations arc obtained from an ex- 
pansion of the thermodynamic fluxes in terms of the cor- 
responding forces. Here, we give only the final equation 
for the relative current between actin and solvent for an 
isotropic active gel [TH] 

ja = -idap- + xdf}(Tlp . (2) 

In this expression ct^^ is are the components of the par- 
tial stress tensor of the gel phase, not to be confused with 
the total stress [18]. It is analogous to the particle-phase 
stress in suspension mechanics [TOll^ In steady state elas- 
tic effects vanish and the stress is purely viscous, so that 
= 2riVg_ai3 — C^M'^ct/S where 77 is the gel viscosity and 
i^g.a/g = 5 {daVg,i3 + dpVg^a) and C^M gives the magnitude 
of the actively generated (isotropic) stress in the gel phase. 
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The first term in expression Q accounts for the diffusive 
current due to gradients in the relative chemical potential. 
The second term describes the fiux that results from gra- 
dients in the gel stress. Finally, the permeation constant 
PgX~^ ^ Vs/^'^j where ^ is the gel's mesh size and r]s the 
solvent viscosity. 

As the only vector for a non polar system is the gra- 
dient vector, within the Onsager linear approach the ac- 
tive contribution to the current is proportional to daAfi 
and vanishes if Afi is constant in space, which we assume. 
Fluxes stemming from the system's activity are captured 
by the last term of the stress expression. The coupling 
parameter ^ here depends on the gel density pg and is 
negative due to the contractility of the motors. 

We proceed using the same approximation as in Ref. [18| : 
The osmotic pressure 11 satisfies the Gibbs-Duhem equa- 
tion dn = Pgdfi, if the volume of the system is indepen- 
dent of the composition and incompressible. We define 
an effective osmotic pressure including active effects as 
77 = 77 + (Ap and, we thus get 

— {Vg,a - Va) = 2l-ldf:sVg^ai3 - i9q7T, (3) 

X 

where we have expressed the gel stress in terms of the 
gel shear rate. There are two competing dissipative mech- 
anisms, the gel viscosity and permeation of the solvent 
through the actin gel. The comparison between these two 
types of dissipation defines a permeation length ~ 
(^Pg ^x)^^^ ^ {''l/'^sY^'^S,- The relative permeation current 
on the left hand side of Eq. ([3]) is negligible compared to 
the viscous dissipation term if we consider the dynamics 
on length scales smaller than the permeation length scale 
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Lp. For typical viscosities rj IQ^rjs and typical mesh 
sizes of ^ of a few tens of nanometers, the characteristic 
length is of the order of a few hundreds of microns and 
thus macroscopic. We therefore neglect permeation in the 
following. 

From now on, we consider only the gel density pg and 
the corresponding velocity field Wg. Dropping the g indices, 
the equations read 

dtp + dapva = -kdP 
2vdpvc.p - dc.n{p) = 0. (4) 

These equations are identical to those of a hydrodynamic 
theory for an effective one component compressible active 
gel HIS]. 

3 Active pre-wetting 

We consider a situation where the gel is assembling at 
the surface located at z = and assume invariance under 
translations in the surface plane. This leaves us with a one- 
dimensional problem. The dynamic equations Q read for 
z > 0: 

dtp + dzpv = -kdP 
r^d.v - n{p) = 0, (5) 

where v denotes the z-componcnt of the gel velocity. The 
second equation follows from integrating Eq. (U). Let us 
recall that the boundary condition on the gel current at 
the surface is pv\z=o = po^p- 

To complete our description, we need to provide an 
expression for the effective pressure 77. Usually, the os- 
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motic pressure iT is a monotonic increasing function of 
the actin density. However, the active contribution is neg- 
ative and has a maximum as a function of density for 
a given crosslink density |21) . The passive osmotic pres- 
sure always dominates at large densities but if activity is 
large enough, the effective osmotic pressure can become a 
non monotonic function of density, which is indicative of 
a phase separation in the solution induced by the contrac- 
tility of the molecular motors. In the following, we use 

n = ap^ + bp^ , (6) 

where the coefficient a depends on activity. For vanishing 
activity, standard three body interactions should domi- 
nate and a should be positive. For large activity, contrac- 
tility should dominate and a should be negative. In prin- 
ciple, the expansion in powers of p should contain linear 
and quadratic terms as well. Whether contractility should 
be refiected in the quadratic or in the cubic term depends 
on how the motor density relates to the actin polymeric 
density. For example, one could argue that contractility 
requires a pair of actin filaments as well as myosin, so if 
the bound myosin concentration is proportional to that of 
the filamentous actin, one would expect a cubic p depen- 
dence. The choice of the functional form is not essential, 
what really matters is the non monotonicity of n{p). If 
a < then the activity leads to contractile stresses. 

We solve Eq. ([5]) in steady state. Combining the two 
equations, we get for the velocity field 

Vy=[^) \fip) , (7) 



leading to 
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where we have introduced the auxiliary function where is a constant. The density thus eventually decays 

exponentially. If po is small, we have On « po and 

f{p) = -karj - ap^ - bp". (8) 

dp kd , ^ 

-T-^ -P , (12) 

The general solution to the steady state equations can az Vp 

then be written in the form such that the characteristic length is Vp/ka- For large val- 

dp 1 9 s ( kdV , ,1 ues of poj the characteristic length is instead VpPo/kdPn 

dz por]Vp [Jp^ P'fip') J 

with 

-z;pexp<j - / -^yjY^j-^dp \ . (10) 

The profile can thus be a simple exponential or it can 
present a "shoulder" presaging the existence of a well- 
defined cortical layer for high contractile activity. In the 
latter case, the profile presents two inflection points, which 
arc determined by d{dzp) / dp — 00. This equation has 
two solutions for —a large enough but still smaller than 
— Oc. In this case, we can define the layer thickness L via 
the condition of mass conservation in steady state p^Vp ^ 
kdPrriL, where pm is approximately given by the value of 
p maximizing /. 

Examples of density profiles for a subcritical activity 
and two different values of po are given in Fig. [5) 

3.2 The large contractile activity regime 



P ' ' \ Jpo p'fip 

From these solutions we infer that the gel density p 
jumps to zero at densities fulfilling f{p) = 0. Indeed, equa- 
tion (fTU|) implies v > with w = for / = and no mate- 
rial is transported beyond such a point. Consequently, as 
soon as / presents zeros a well-defined cortical layer with 
a sharp edge can be formed. The existence of zeros of / 
depends on the activity. There exists a critical value Oc 
beyond which a wetting layer appears. For —a > — > 
the function has two zeros, while for —a < — Oc it has 
none. If the activity is equal to the critical value Oc, the 
function / has exactly one zero. As we will discuss be- 
low, the transition at a ~ ac can be viewed as an active 
or non-equilibrium analog of a pre- wetting transition |14l 

Eg. 



We consider now the case —a > —Ac, such that f{p) has 
3.1 The weak contractile activity regime two real roots pi > p2. From Equation © and since v > 



The weak contractile activity regime is defined by —a < 
— flc, such that /(p) < for all p. Since p tends to zero as 
z — 00, we have in this case /(p) ~ —kdT] and thus 



we infer that for a density p{z = 0) = po < P2 the 
density approaches with increasing z. Indeed, in this case 
dp/dz < for all z > 0. As a consequence, the density 
profile behaves very similarly to the weak activity case 



^__J^2^^ /-In— 1 (11) 

dz PoVp^ y ^ Po J ' ^ ^o^c- In contrast, for po > p2, the density approaches 
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a) p 



b) V 




Fig. 2. Steady state density and velocity profiles in the weakly 
active regime, —a < —ac from numerical solutions of Eqs. ©. 
a) Densities as a function of the distance from the polymerizing 
surface for surface densities po = 1.3 and po = 2.0, respectively. 
Inset: Auxiliary function /(p). b) Corresponding velocity pro- 
files; dotted: po ~ 1.3 and solid: po ~ 2.0. Parameters are 
a = —karj and b = Q.hkdrj. 

the value pi. As w — > for p — > pi, this value is reached 
at a finite distance z from the surface. 

We can calculate the density profile for p « pi. Then 



Hp) 



-a{p — pi) with a = |/'(pi)| such that 



pa 



P'fip') 



kdV \P- Pi\ 
\Po - Pi I 



api 



(14) 



and 



dp ap\ 
dz por]Vp 



iP~ Pi) 



\P- Pi\ 



1 -k^rj/api 



\Pa - Pi I 



(15) 



We can solve the latter equation explicitly. In the case 
Po ~ pi we get for the whole density profile 



p(z) = pi + (po - pi) <{ — — (L 
. Po Vp 



api/kdV 



(16) 



Note, that the density profile in the vicinity of L depends 
on the value of api/kdJj that fixes the slope of the profile 
when p = pi . In Figures [3] and IH we present steady state 
solutions for various nucleator densities po for api/kdf] < 
1 and api/kdTj > 1, respectively. 

From the approximate density profile (jl6p , we estimate 
the cortex thickness L by setting z = 0, yielding 



L 



PQ_Vp_ 

Pi kd 



(17) 



As can be seen in Figs. [3] and SI for p w po it approximates 
well the thickness obtained from the numerical solution to 
Eqs. ([5]). For po > pi the length increases linearly with po, 
however, with a slope that is different from the one given 
in Eq. (d?]). 

4 Discussion 

In this work, wc have introduced a hydrodynamic theory 
for describing the dynamics of active gels in presence of fil- 
ament polymerization and depolymerization. Our descrip- 
tion of filament assembly and disassembly requires some 
comments. Specifically, we have taken the depolymeriza- 
tion rate and viscosity to be constant. Our main results are 
more general and do not depend on this assumption. The 
same behavior is obtained if we keep the gel density depen- 
dences of depolymerization rate and viscosity. In general, 
however, the depolymerization process is more complex as 
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a) p 



Fig. 3. Steady state density and velocity profiles in the 
strongly active regime, —a > — Oc, from numerical solutions 
of Eqs. (O. a) Densities as a function of the distance from the 
polymerizing surface for surface densities po = 1.2, po ~ 1-3, 
and Po = 2.3, respectively. Inset: Auxiliary function /(p), 
which has zeros at pi = 2.08 and p2 = 1.23. b) Corresponding 
velocity profiles; dotted: po — 1.2, solid: po — 1.3, and dashed: 
po ~ 2.3. Inset: cortex width as a function of po from numerical 
solution of Eqs. ((5|) (dots) and from Eq. (|17p (line). Parameters 
are a = —1.15kdrj and b = O.Skdrj. 

the rate of filament disassembly depends on mechanical 
stresses, the degree of severing, and the presence of actin 
associated proteins like cofilin or gel-solin. In addition the 
fact that filaments depolymerize from their ends might in- 
troduce a gradient in the effective depolymerization rate. 
Together these processes might even lead to filament sub- 
populations virith different turnover rates [23] • Similarly, 
the polymerization process is more involved than assumed 




Fig. 4. Steady state density and velocity profiles in the 
strongly active regime, —a > — Oc, from numerical solutions 
of Eqs. ((Sjl. a) Densities as a function of the distance from the 
polymerizing surface for surface densities po = 0.9, po ~ 0.94, 
and Po = 1.4, respectively. Inset: Auxiliary function f{p), 
which has zeros at pi = 1.0 and p2 — 0.93. b) Correspond- 
ing velocity profiles; dotted: po = 0.9, solid: po — 0.94, and 
dashed: po = 1.4. Inset: cortex width as a function of po from 
numerical solution of Eqs. ([5|) (dots) and from Eq. (|17|) (line). 
Parameters are a = —A.bkdr] and b — S.Skdrj. 

in the present work. We have restricted it to be spatially 
localized at the membrane, while actin filaments also grow 
in the bulk and can be nucleated away from the mem- 
brane. This possibility will be discussed in a forthcom- 
ing publication. Wc believe, however, that linking cortex 
formation to an out of equilibrium wetting phenomenon 
opens a new way of thinking about the problem worth 
being fully investigated. This picture is a non-equilibrium 
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